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Abstract 


This  report  presents  detailed  derivations  of  the 
expressions  for  the  steady-state  disturbance  velocity  at  a  pro¬ 
peller  blade  due  to  pressure  loading  and  thickness.  A  procedure 
for  separating  lifting  line  velocities  from  the  total  is  outlined, 
and  conclusions  are  drawn  regarding  the  presence  of  camber  and 
incidence  corrections  for  propellers  with  symmetrical  blades  and 
chordwise  load  distributions. 
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1.  Introduction 


A  linearized  theory  for  propellers  can  be  developed  as  a  logical 
extension  of  the  theory  of  finite  wings  in  incompressible  flow.  The 
principal  difference  is  that  a  point  on  a  propeller  blade  describes  a 
helicoidal  path  as  it  moves  through  the  fluid,  while  a  point  on  a  wing 
simply  moves  along  a  straight  line.  Although  the  general  form  of  the 
expressions  relating  shape  to  pressure  distribution  is  similar  in  both 
cases,  the  propeller  expressions  are  naturally  more  complicated  due  to 
the  more  involved  geometry  of  the  motion. 

We  begin  by  stating  the  usual 'assumptions  that  the  fluid  is 
frictionless,  incompressible,  free  of  cavitation  and  infinite  in  extent. 
The  propeller  is  considered  to  be  operating  in  an  axially  directed  stream 
whose  magnitude,  is  a  function  of  radius  only,  and  it  is  assumed  that  there 
are  no  extraneous  solid  boundaries.  The  flow  is  therefore  steady  relative 
to  a  coordinate  system  rotating  with  the  propeller. 

The  propeller  is  assumed  to  have  K  symmetrically  spaced  iden¬ 
tical  blades.  In  the  present  work  the  existence  of  a  propeller  hub  is 
neglected  entirely.  The  propeller  blades  begin  at  some  radius  r^  which 
corresponds  to  the  radius  of  the  hub  in  the  actual  propeller.  Both  the 
inner  and  outer  extremities  of  the  blades  are  therefore  regarded  as  free 
ends  when  the  lift  is  required  to  be  zero..  This  assua^tion  is  obviously 
not  particularly  realistic,  and  work  is  presently  in  progress  to  represent 
the  hub  properly  as  a  solid  cylindrical  boundary.  However,  since  the 
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inner  part  of  the  blade  near  the  hub  contributes  only  a  small  portion 
of  the  total  propeller  thrust,  the  "hubless  propeller"  approximation 
may  be  adequate. 

In  the  linearized  case  we  assume  that  the  disturbance  to  the 
flow  caused  by  the  propeller  is  small,  which  requires  that  the  blades 
be  thin  and  that  the  camber  and  incidence  of  the  blade  sections  relative 
to  the  undistrubed  flow  be  small.  As  in  linearized  wing  theory,  the 
effects  of  thickness  and  of  camber  and  incidence  can  be  considered  sep¬ 
arately.  Both  can  be  represented  by  singularities  distributed  on  stream 
surfaces  formed  by  the  undisturbed  relative  flow  past  the  blades. 

Blade  thickness  can  be  represented  by  a  distribution  of  sources 
and  sinks  whose  strength  is  proportional  to  the  slope  of  the  thickness 
form  in  the  streamwise  direction.  Similarly,  the  discontinuity  in  pressure 
across  the  lifting  surface  can  be  generated  by  a  distribution  of  vor- 
ticity  whose  axis  is  on  the  lifting  surface  and  normal  to  the  oncoming 
flow.  These  are  designated  "bound  vortices."  The  condition  of  contin¬ 
uity  of  vorticity  requires  the  presence  of  an  additional  distribution  of 
vorticity  whose  axis  is  oriented  in  the  direction  of  the  relative  flow. 
These  vortices  known  as  "trailers"  are  present  not  only  on  the  blade 
surface  but  extend  to  infinity  in  the  wake  behind  each  blade. 

There  are  obviously  a  large  number  of  possible  types  of  problems 
to  which  steady-state  propeller  theory  may  be  applied.  However,  the 
following  three  are  of  principal  interest: 

1)  Determine  the  shape  of  the  blade  sections  for  a  prescribed 
radial  and  chordwise  load  distribution. 

2)  Determine  the  maximum  camber  and  incidence  of  a  given  type 
of  section  at  each  radius  in  such  a  way  that  a  prescribed 
radial  load  is  achieved  with  the  sections  operating  at  their 
ideal  angle  of  attack. 

3)  Determine  the  load  distribution  for  a  given  shape. 
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The  first  problem  is  the  most  straight  forward.  We  assume  that  the  blade 
outline  and  thickness  distribution  have  been  determined  from  considerations 
of  strength  and  cavitation.  With  both  the  load  andthickness  specified, 
the  strengths  of  the  sources,  bound  vortices  and  trailers  are  immediately- 
known.  The  velocity  normal  to  the  blades  at  each  point  induced  by  each 
of  these  three  singularity  distributions  may  be  determined  by  integration. 
The  slope  of  the  section  mean  line,  according  to  the  linear  theory,  is 
equal  to  the  normal  component  of  the  disturbance  velocity  divided  by  the 
magnitude  of  the  undisturbed  approach  flow.  Consequently,  if  the  normal 
component  of  the  disturbance  is  known,  the  section  shape  can  easily  be 
determined  by  integration. 

An  important  consideration  is  that  the  shape  of  the  section 
mean  lines  required  for  a  particular  load  depends  not  only  on  the  load 
but  also  on  the  thickness  of  the  blades.  This  is  not  true  in  the  case 
of  a  planar  lifting  surface  since  the  sources  cannot  induce  a  velocity- 
normal  to  the  plane  on  which  they  are  situated.  However,  in  the  case 
of  a  propeller,  the  sources  can  induce  a  velocity  normal  to  the  blade 
on  which  they  are  situated  as  well  as  normal  to  the  adjacent  blades. 

The  second  problem  is  somewhat  more  complicated  since the  chord- 
wise  xi ad  is  initially  unknown.  However,  as  in  the  case  of  wing  lifting 
surface  theory,  the  chordwise  load  can  be  assumed  to  be  composed  of  a 
finite  number  of  modes  whose  form  is  determined  from  two -dimensional 
theory.  The  normal  velocity  induced  by  each  mode  can  be  determined  in 
the  same  way  as  in  the  first  problem.  Since  the  total  load  at  each 
chordwise  section  is  prescribed,  the  amplitudes  of  some  of  the  chordwise 
modes  are  known  initially.  The  camber,  incidence  and  amplitudes  of  the 
remaining  modes  can  then  be  obtained  by  collocation. 
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The  third  case,  which  is  sometimes  call  the  "inverse  problem," 
can  be  solved  by  assuming  that  the  unknown  load  distribution  is  composed 
of  a  double  summation  of  chordwise  and  radial  modes.  The  normal  velo¬ 
city  induced  at  a  set  of  points  by  each  mode  can  be  determined  as  before 
and  the  unknown  amplitudes  obtained  by  collocation.  This  type  of  problem 
would  arise  if  one  wished  to  analyze  the  performance  of  a  propeller 
operating  at  other  than  design  conditions. 

In  the  following  sections,  detailed  derivations  are  given  for 
the  disturbance  velocities  due  to  the  source  and  vortex  systems  repre¬ 
senting  the  propeller.  Much  of  this  material  may  be  found  in  recent 
publications  by  a  number  of  authors.  In  particular,  it  should  be 
mentioned  that  the  results  given  here  in  section  4,  5/  and  6  are  in 
complete  agreement  with  corresponding  results  given  by  Sparenberg  [  l] 
and  by  Plan  12],  provided  the  proper  conversion  of  nomenclature  is 
made.  Numerical  techniques  and  results  have  not  been  included  in  this 
report  but  may  be  found  in  references  [3]  -  [7] . 
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2.  Geometrical  Considerations 


We  begin  by  defining  a  Cartesian  (x,  y,  2)  coordinate  system 
which  is  fixed  on  the  propeller.  As  shown  in  Fig.  2.1,  the  x  axis 
is  the  axis  of  revolution  with  positive  distances  measured  downstream. 

The  y  axis  is  selected  so  as  to  pass  through  the  tip  of  one  blade,  while 
the  z  axis  completes  the  right-handed  system.  It  is  also  convenient  to 
define  a  cylindrical  (x,  r,  9)  system  where  the  x  axis,  as  before,  is  the 
axis  of  revolution.  The  radial  coordinate  is  denoted  by  r,  and  the 
angular  coordinate  (measured  clockwise  starting  from  the  y  axis)  by  0. 

The  equations  relating  the  two  systems  are: 

_  r 2  2 

y  =  r  cos 9  r  =  Jy  +  z 

(2.1) 

z  =  r  sinG  0  =  tan-1(z/y) 

In  order  to  relate  corresponding  points  on  each  of  the  K 
blades,  we  define  6^  as  the  G  coordinate  of  the  point  at  the  tip  of 
the  k'th  blade.  Assuming  that  the  blades  are  symmetrically  arranged, 
these  angles  are: 

6  =  k  =  1,  2...K  (2.2) 

K 

which  is  illustrated  in  Fig.  2.2. 

In  order  to  conserve  symbols,  we  will  assume  at  the  outset 
that  the  propeller  radius  R  is  equal  to  unity.  This  is  equivalent  to 
saying  that  all  length  dimensions  (including  those  present  in  linear 
velocities)  have  been  divided  by  R. 

The  undisturbed  flow  past  the  blades  is  composed  of  a  tan¬ 
gential  component  car  due  to  the  rotation  of  the  propeller,  and  an 
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VIEW 


FIG. 2.1  COORDINATE  SYSTEMS 


axial  component  VA(r)  due  to  the  oncoming  stream.  Since  there  is  no 
radial  component,  the  streamlines  lie  on  cylindrical  surfaces,  r  =  con¬ 
stant.  The  angle  between  the  undisturbed  relative  flow  and  the  yz  plane 
can  be  seen  from  Fig.  2.3  to  be: 

3(r)  =  tan-1  (-£ - )  =  tan_1(^I)  (2.3) 

sur 

where  g  is  known  as  the  advance  angle  and 


X(r)  =  r  tang(r) 


(2.4) 


is  defined  as  the  advance  coefficient.  The  streamlines  are  therefore 
helices  whose  pitch  P(r)  is  equal  to  2n)(r). 

The  sine  and  cosine  of  the  advance  angle,  which  will  appear 
frequently  in  the  following  sections,  can  be  expressed  in  terms  of  r 
and  X 

sinp(r)  -  ^ 

V?+  ^(r) 

(2.5) 


coag(r)  -  ~ —  1 ' 

+  ^(r) 

as  can  easily  be  seen  from  (2.3). 

We  next  assume  that  the  propeller  blades  lie  approximately  on 
the  helicoidal  surfaces  swept  out  by  the  undisturbed  flow  past  the 
radial  lines 


x  -  0  j 

This  does  not  impose  any  restriction  on  the  blade  outline,  but 
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neglects  the  effect  of  blade  rake. 

The  k'th  blade  and  Its  wake  Is  therefore  located  approximately 
on  the  helicoidal  surface 

H^x,  r,  0)  -  x  -  X(r)(e  -  6k)  -  0  (2.7) 

It  is  also  convenient  to  define  an  orthogonal  curvilinear 
coordinate  system  (s,  n,  r)  on  each  of  the  surfaces.  The  s  coordinate 
is  formed  by  the  intersection  of  an  axial  cylinder  and  the  surface  H^, 
and  Is  therefore  tangent  to  the  streamlines  of  the  undisturbed  flow.  The 
r  coordinate  is  radial,  as  before,  and  the  n  coordinate  is  perpendicular 
to  r  and  s  is  directed  in  such  a  way  that  it  has  a  positive  axial 
component.  This  notation  is  illustrated  in  Pig.  2.1. 

If  the  pitch  of  the  helicoidal  surface  is  independent  of  radius, 
the  r  coordinate  through  any  point  on  remains  on 
is  therefore  normal  to  Hj^.  If  the  pitch  is  a  function  of  radius,  this 
is  not  necessarily  true  so  that  n  may  depart  slightly  from  the  true  nor¬ 
mal  to  the  helicoidal  surface.  However,  it  is  assumed  that  variations  in 
pitch  are  sufficiently  gradual  for  this  discrepancy  to  be  negligible. 

As  can  be  seen  from  Pig.  2.1,  the  expressions  relating  the 
(s,  n,  r)  and  (x,  r,  9)  coordinate  systems  are: 

X(r)  x  +  r2(0  -  6  ) 
s  «  x  sing  +  r(8  -  6.)  cosg  ■  — - - —  ■  ■  ■■■  - 

Vr2.  X2(r) 

(2.8) 

xr  -  rX(e  -  6  ) 

n  =  x  cosg  -  r(0  -  6.  )  sing  -  -  —  . — — 

Vr2  ♦ 


H.  .  The  n  coordinate 
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On  any  one  of  the  K  helicoidal  surfaces,  the  relationship  between  x  and 
6  can  be  obtained  by  combining  (2.7)  and  (2.8) 


s  =  >Jr2  +  X2(r)  (0  -  4k)  (2.9) 

We  can  now  define  the  blade  outline  and  the  shape  of  the  blade 

sections  in  the  (s,  n,  r)  system.  The  s  coordinates  of  the  leading  and 

trailing  edges  are  s  (r)  and  s  (r)  respectively,  while  the  angular 
L  1 

coordinates  are  designated  ^(r)  and 

=  sT(r)  -  sL(r)  (2.10) 

The  blade  sections  can  be  decomposed  (as  in  two-dimensional 
wing  theory)  into  a  symmetrical  thickness  form  t(s,  r),  a  mean  line 
c(s,  r)  and  an  angle  of  incidence  o(r).  The  maximum  thickness  of  a 
section  at  a  particular  radius  is  denoted  tQ(r)  and  the  maximum  camber 
is  cQ(r).  This  notation  is  shown  in  Fig.  2.4. 
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sL(r) 
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3.  normal 

According  to  linear  theory,  the  thickness  of  the  blades  can  be 
generated  by  a  distribution  of  sources  and  sinks  located  on  the  a  col  dal 

surfaces  representing  the  blades.  We  consider  first  a  point  source  of 
unit  strength  located  at  a  point  with  Cartesian  coordinates  (5,  T],  Q). 

The  velocity  potential  if  this  source  Is 


y>  z)  -  — -  1  £=  - .  ■  -7  (3.1) 

4tY(x  -  if*  (y  -  ti)2+  (z  -  c f 
and  the  Cartesian  components  of  velocity  at  (x,  y,  z)  are 


u(x,  y,  z,  g,  TV,  C)  =  M  . . 

4r[(x  -  g)2+  (y  -  H)2+  (z  -  C)2]3/2 

(3.2) 

v(x,  y,  z,  5,  n,  c)  =  U  - - - 

4*  [(x  -  g)2+  (y  -  T))2+  (z  -  02J3/2 


w(x,  y,  Z,  i,  n,  c)  -  |£  - - £-=-i - 

4"  [(x  -  g)2+  (y  -  n)2+  (z  -  O2]3/2 

The  bars  on  the  symbols  (u,  v,  w)  denote  velocities  due  to 
a  point  source  of  unit  strength.  The  total  velocity  at  a  point  due  to 
a  distribution  of  sources  is  identified  by  the  same  symbol  but  without 
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a  bar  on  top*  It  should  be  mentioned  that  the  symbols  (u,  v,  w) 
will  not  be  used  exclusively  for  source  velocities  but  will  be  used  to 
denote  disturbance  velocity  components  In  general.  Since  we  are  consid¬ 
ering  only  sources  in  this  section,  there  is  no  need  for  further 
identification.  However,  additional  symbols  will  have  to  be  introduced 
later  to  distinguish  between  disturbance  velocities  due  to  source  and 
vortex  distributions. 

The  velocity  components  given  in  (3*2)  can  be  converted  to 
cylindrical  coordinates  using  the  equivalent  of  (2.1 ) 

y  =  r  cosS  z  =  r  sine 

„  (3.3) 

T|  =  p  coscp  C  =  P  sincp 

where  (x,  r,  0)  is  the  point  where  the  velocity  is  to  be  evaluated  and 
(?>  P>  <p)  is  the  location  of  the  source.  The  axial  tangential  and 
radial  components  are 

ua(x,r,0,5,p,!p)  - - - - 

4n£(x  -  ?)2+  r2+  p2-  2rp  cos(tp  -  0)J^2 

ut(x,r,0,5,p,cp)  =  w  cos©  -  v  sin0  - - sin(cp  -  9) _ 

^[(x  -  §)2+  r2+  P2-  2rp  cos(cp-0)J^/2 

ur(x,r,9,5,p,q>)  =  »  sine  +  v  cos0  =  - r  -  pcos(cp  -  e)  _ 

4n^(x  -  §)2  +  r2  +  p2-  2rpcos(<p-0)J^/2 


lb 


The  velocity  components  given  in  (3**0  can  be  resolved  into 
components  in  the  helicoldal  coordinate  system  by  means  of  (2.8) 


us(x,  r,  9,  §,  p,  9) 


^(x,  r,  0,  5,  p,  cp) 


X(r)y  r5t 

>/r2  +  A2(r) 
"/r2  +  X2(r) 


(3.5) 


In  this  application,  we  are  only  interested  in  the  velocity 
normal  to  the  helicoidal  surface.  This  can  be  obtained  by  combining  (3«M 
and  (3.5) 


{^(x,  r,  9,  §,  p,  <p) 


r(x  -  g)  +  A(r)p  sin(cp  -  9) 

WVr2  +  A2(r)£(x  -  ?)2+  r2+  p2-  2rp  cos(cp-0)p/2 


(3.6) 

For  points  on  the  helicoidal  surfaces  representing  the  blades, 
x  and  5  are  related  to  9  and  9  by  (2.7).  If  we  adopt  the  point  of  view 
that  the  velocity  is  always  to  be  computed  at  a  point  (x,  r,  0)  on  the 
first  blade,  while  the  position  of  the  source  is  at  some  point 
(5,  p,  9+  6^)  on  the  k'th  blade,  there  follows 


x  =  X(r)0 
§  =  X(p)9 


(3.7) 
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so  ifaftt  (3 .6)  becomes 


r(A(r)6  -  X(p>p)  +X(r)p  ain(9  +  0 


\(r,  0,  p,  9,  k) 


inWr2+ A2(r)[(*(r)0  -  *(p>p)2+  r2+  p2-  2rpoos(^+«k-| 


(3.9) 

If  X  is  independent  of  radius ,  the  term  (X(r)0  -  A(p)«p)  appearing  both 
in  the  numerator  and  denominator  of  (3*8)  can  be  replaced  hyX(0  "  ?)» 

In  this  case,  the  angles  9  and  0  in  (3.8)  appear  only  in  the  combination 
(<p  -  e),  which  is  equivalent  to  saying  that  the  velocity  depends  only  on 
the  angle  between  the  source  and  the  point •  Since  this  represents  a 
great  simplification,  it  is  worthwhile  in  the  general  case  to  introduce 
the  approximation 

A(r)e  -  X(p)<p  -  X(p)(<p  -  e)  =  -  A(p)u  (3*9) 

where,  by  definition 

H  *  cp  -  0  (3*10) 

»Ph«.  geometrical  interpretation  of  this  approximation  is  that  the 
axial  distance  between  the  point  at  (x,r,e)  and  the  singularity  at  (§,P,<p) 
is  assumed  to  be  the  same  as  if  the  pitch  of  the  helicoidal  surface  at  p 
were  the  same  as  at  r.  If  the  variation  of  X  with  r  is  gradual,  the  error 
in  distance  will  be  snail  when  p  is  close  to  r.  Consequently,  the 

error  introduced  in  determining  the  velocity  induced  by  nearby  source 
elements  would  presumably  be  small.  For  distant  elements,  on  the  other 
hand,  a  change  in  the  axial  component  of  the  distance  between  two  ele¬ 
ments  will  not  have  a  large  effect  on  the  direction  or  magnitude  of  the 
Induced  velocity. 
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Of  course,  when  X  Is  constant,  there  is  no  approximation  involved.  ’’ 
Substituting  (3*9)  into  (3.8)  and  sensing  orer  the  K 'blades, 
ve  obtain  the  result 

S(r,p,n)  - - -AlP,).-,.  ^  - l3l  +,J>  sl?k  - 

lm'/r2+  X2(r)  2^  [X2(p)^2+  r2  +  p2  -  2rp  cosfo  +  C^Jp/2 

(3.11) 

where  S(r,p,p,)  represents  the  velocity  normal  to  the  helicoidal  surface  at 
a  point  (r,  6)  on  the  first  blade  induced  by  K  unit  point  sources  at  a 
radius  p.  The  angular  coordinates  of  the  sources  are  (cp  +  -  8) 

measured  from  the  y  axis,  or  +  8^)  measured  from  the  point  (r,  8). 

It  can  be  seen  from  (3.11)  that  S(r,p,u)  is  an  odd  function  of  n, 

S(r,p,-ji)  =  -S(r,p,(i)  (3.12) 

This  is  obviously  true  when  k  *  1.  When  k  >  1,  the  contribution  of  the 
k'th  blade  to  S  is  odd  provided  that  we  change  the  sign  of  5^  as  well 
as  u.  Since  the  sign  change  is  equivalent  to  summing  the  blades  in 
the  reverse  order,  the  total  result  is  unaffected,  hence  ve  conclude  that 
S  is  odd. 

We  next  need  to  determine  the  strength  of  the  source  distribution. 
According  to  linear  theory,  the  source  strength  per  unit  area,  o,  is 
given  by^ 


o(p,  tp)  -  VQ(p)  H  (p,  qp) 


(3.13) 
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\i. 

where  VQ  Is  the  magnitude  of  the  undisturbed  relative  flow  and 

Is  the  derivative  of  the  thickness  form  in  the  streamwlse  direction. 
It  is  evident  from  Pig.  2.3  that 


y»>  _  yp)/pa*  xg(p) 

sinp(p)  "  ^(p)  (3*l4) 


so  that 


VA(P)  'Ip2  +  A2(p)  §§(p,<p) 

A(p) 


(3.15) 


The  total  normal  velocity  at  any  point  on  the  first  blade  can  now  be 
obtained  by  combining  (3«ll)  and  (3*15)  and  integrating  over  the  sur¬ 
face  of  the  first  blade 

1.0  sJp) 

(r,e)  =  J  J  o(p,cp)  S(r,  p,  n)  ds  dp  (3-l6) 

p=rh  s=sL(p) 


The  chordwise  integration  in  (3*16)  can  be  expressed  entirely  in  terms 
of  the  variable  u  by  introducing 


ds 


tp  =  |i  +  e 


(3.17) 
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so  that 


While  the  Integration  In  (3«l8)  cannot  he  carried  out  explicitly  in 
most  cases ,  various  numerical  schemes  are  available  which  will  yield 
results  of  sufficient  accuracy.  It  is  Important  to  note  that  the 
Integral  In  this  case  is  not  singular  since  S(r,p,n)  has  a  finite 
limit  when  r-*p  and  4-0. 
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4.  Distribution  of  Bound  Vortices  and  Trailers 


The  pressure  loading  on  the  blades  can  be  represented  by  a 
distribution  of  radially  oriented  bound  vortices  lying  on  the  helicoldal 
surfaces  representing  the  blades.  We  again  Introduce  the  duany  cylindrical 
coordinates  (S,P,cp)  to  denote  the  location  of  an  element  of  the  vortex 
sheet,  while  (x,r,  0)  are  the  coordinates  of  a  point  at  which  the  velocity 
is  to  be  determined.  The  strength  of  the  bound  vortex  sheet  per  unit 
of  length  along  the  helix  at  radius  p  is  y(  P,«p).  The  pressure  difference 

/Q\ 

across  the  surface  according  to  linear  theory  isv  ' 


Ap  (p,q>)  =  P  V0(p)  Y(P*<P)  (4.1) 

where  p  Is  the  fluid  density.  Substituting  (3.14)  for  the  approach 
velocity  Vq  In  (4.1 )  there  follows 

,  v  pva(p)/p2  +  X2(p)  y(p,«p) 

^P(P><P)  - - 

*(p)  (4.2) 

The  total  circulation  around  a  blade  at  radius  p  Is 
sT(p)  «T(p) 

T(p)  =  \  Y(p»<P)ds  - 

SL(P) 

where  0^  and  8p  are  the  angular  coordinates  of  the  leading  and 
trailing  edges,  respectively,  and  the  element  of  arc  length  ds  has  been 


\  y(Pj'P)  J" 

VP) 


*2(p)cUp 


(4.3) 
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expressed  In  terns  of  dtp  using  (2.9) 


In  order  to  preserve  continuity  of  vortlcity,  there  mist  also 
he  a  system  of  trailing  vortices  whose  axis  la  In  the  a  direction.  The 
requirement  of  continuity  of  vortlcity  on  a  hellcoidal  surface  can  be 
obtained  very  easily  by  keeping  track  of  the  total  vortlcity  entering  a 
differential  element  of  the  surface,  as  shown  In  Fig.  4.1.  The  strength 
of  the  vortex  lines  entering  the  four  sides  of  the  element  Is 


7  2  g 

p  +  (p)dtp 

bound  circulation  entering  top  =  -Y(P+dp>cp)  \f  (p+dp)2+  A2(p+dp)d«p 
trailers  entering  left  =  Y_(p,9)  dp 

(4.4) 

trailers  entering  right  -  -Y*(p><p  +  dtp)  dp 
Setting  the  sum  equal  to  zero  gives  the  result 


(4.5) 


where  ys  denotes  the  strength  of  the  trailer. 

Equation  (4.4)  holds  only  for  points  In  the  Interior  of  the 
blades.  At  the  leading  and  trailing  edges,  a  trailer  Is  generated,  so 
to  speak,  whenever  a  bound  vortex  tries  to  run  off  the  edge.  This  can 
he  seen  by  taking  as  the  control  area  an  infinitesimal  triangular  element 
as  shown  in  Fig.  4.2.  In  order  for  the  blade  edge  to  form  the  hypotenuse 
of  the  triangular  element,  the  relationship  between  dxp  and  dp  must  be 


8.  -  JJ-4P 


at  the  leading  edge. 

(*.6) 


“Sp“d<> 
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at  the  trailing  edge 


FIG.  4.1  ELEMENT  OF  VORTEX  SHEET 
IN  INTERIOR  OF  BLADE 


At  the  leading  edge,  the  bound  circulation  entering  the  bottom  as  before 
is 

y  (p,^)/p2+  >?(P)*P  (MO 

while  in  this  case,  there  is  none  leaving  the  top.  This  must  be  balanced 
by  the  total  strength  of  the  trailers  entering  the  right  side,  which  to  a 
first  order  is 


-YS(P^)  dp 

Combining  this  with  (4.6)  and  (4.7)  gives  the  result 


(4.8) 


ys(p,Ql)  =  yCp,^) 


(4.9) 


A  similar  analysis  of  circulation  entering  a  triangular  element  at  the 
trailing  edge  shows  that  Y  must  also  be  discontinuous  at  that  point. 

If  Ys(p,0^)  is  the  strength  of  the  trailer  Just  inside  the  trailing 
edge,  and  Ys(p>g£)  is  the  strength  just  outside  the  trailing  edge,  the 
following  expressions  hold  to  a  first  order 

j  g  2 

bound  circulation  entering  bottom  =  y(p>9q,)  v  P  +  X  (p)  dtp 


bound  circulation  leaving  top  =  0 

(4.10) 

trailers  entering  left  =  Yb(P>6l)  dp 
trailers  entering  right  =  -Ys(P>fi£)  dp 
Setting  the  total  equal  to  zero  and  introducing  (4.6),  there  follows 


y8(p,^)  -  Ys(p,ef) 


-y(P/$ji)  \f 


2 

P  + 


A2(p) 


ST 


(4.U) 
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FIG.  4.2  ELEMENT  OF  VORTEX  SHEET 
AT  TRAILING  EDGE 


Bm  total  strength  of  the  trailers  st  any  point  is  obtained  by 
integrating  (4.$)  and  adding  the  contributions  of  the  edges  froa 
M)  and  (4.11) 


for  <  9 

Noting  that  the  limits  of  integration  above  are  functions  of  p,  it  is 
evident  that  (4.13)  is  equivalent  to 

^(P)  _ 

y,(p,<p)  ■  -  4-  J  V(p,<p)  JeS  *  A2(p)  *»>  -  -  ^  (*.H) 


for  6^,  <  9 


25 


5.  Normal  Velocity  Induced  by  Bound  Vortices 

We  determine  next  the  velocity  induced  at  a  point  (x,  r,  6) 
by  a  bound  vortex  of  unit  strength  located  at  a  point  (5,  p,  cp).  Since 
there  is  no  such  thing  as  a  "point  vortex/'  what  ve  will  determine  is 
the  velocity  Induced  per  unit  length  of  the  vortex.  This  can  be  determined 
from  the  law  of  Biot-Savart,  which  in  this  case  can  be  written: 

(u,  v,  v)  -  2L2L2,  (5.1) 

4tt|d|j 

where  r  is  a  unit  vector  in  the  radial  direction,  D  is  a  vector  from 
the  vortex  (5,  p,  cp)  to  the  point  (x,  r,  0),  and  (u,  v,  w)  are  the 
Cartesian  components  of  velocity.  The  components  of  D  can  be  seen 
from  Fig.  5«1  to  be: 

D  .[(a  -  ?),  r  cos0  -  p  coscp,  r  sinB  -  p  sin  cpj  (5.2) 

while  the  components  of  a  unit  vector  in  the  radial  direction  are 

r  =  £0,  coscp,  sincpj  (5.3) 

Substituting  (5.2)  and  (5*3)  in  (5.1)  yields  the  following  expression 
for  the  velocity  components 


u  = 


r  sin(cp  -  fl) 


kn 


[(x  -  §)2  +  r2  +  p2  -  2rpcos(cp  -  e)j 


3/2 
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Equation  (5*6)  expresses  the  Telocity  Induced  at  (x,  r,  |)  by  a  unit 

bound  vortex  at  (?,  p,  cp).  If  ve  again  consider  the  special  case  of 

the  velocity  Induced  at  (r,  9)  on  the  first  blade  by  a  unit  bound  vortex 

at  (p,  cp  +  6.  )  on  the  k'th  blade,  (5.6)  becomes 

k  (5.T) 

-  -r2  sin(cp-0+6j  -  X(r)fx(p)<p  -  A(r)e|cos(«p-t^k) 

n(r,6,P,<p»k)  - t=== r~ - - -■  — — -  -  _ 

kn)/r2+  >?(r)  j^(A(p)cp  -A(r)9)2+  r2+  p2-  2rpcos(cp-0+fik)j 

where  x  and  ?  have  been  eliminated  using  (3* 7)* 

Introducing  the  approximation 


X(r)9  -  A(pX>*  -X(p)  (cp  -  6)  =  -X(p)h  (5.8) 

as  before  and  sunning  over  the  K  blades,  we  obtain  the  result 


B(r,p,u) 


4tt 


1  r2+  X2(r) 


-r2  sln(u  +  6k)  -  A(r)X(p)ucos(|i  +  4k) 
|a2(p)u2  +  r2  +  p2  -  2rpcos(n  +  «k)j3^2 


In  equation  (5«9)>  B(r,  p,  u)  represents  the  normal  velocity  Induced 
by  K  unit  bound  vortices  and  is  analogous  to  the  expression  for  the 
velocity  Induced  by  the  unit  sources  given  in  (3-ll).  It  is  also 
evident  that  B  is  an  cdd  function  of  u  provided  that  the  blades  are 
symmetrically  arranged. 

The  total  velocity  at  (r,  0)  induced  by  a  distribution  of 
bound  vortices  can  now  be  written  as  an  Integral  over  the  surface  of  the 


first  blade 


v(p,<p)  B(r,p,n)  ds  dp 


(*>e) 


■  f 


®m(p) 

Jn 

a-^Cp) 


(5.10) 


Y(p^i  +  e)  B(r,p,n)\/f£  +  >f(p)  dp 


In  this  case  the  Integral  is  singular  at  the  point  (p  =  r,  4  =  0) 
so  that  the  Cauchy  principal  value  oust  be  'taken. 
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6.  Borna!  Velocity  lhduced  by  Trailers 

We  begin  by  determining  the  velocity  induced  at  (x,  r,  |)  by 
a  trailer  of  unit  strength  originating  at  a  point  (§,  p,  qp)  and  extending 
to  infinity  downstream,  in  order  to  distinguish  a  general  point  on 
the  trailer  from  its  starting  point,  the  former  is  identified  by  the 

I  f  I 

dumny  coordinates  (§  ,  p  ,  =P  ).  Since  the  radius  of  the  trailer  is 
assumed  to  be  constant,  the  prime  on  p  is  unnecessary  and  will  therefore 
be  omitted. 

The  velocity  can  be  expressed  according  to  Biot-Savart '  s  lav 
and  an  integral  along  the  vortex 


(u,  v,  w) 


ds ' 


(6.1) 


where  s  is  a  unit  vector  tangent  to  the  trailer  and  ds*  is  a 
differential  element  of  are  length  along  the  trailer.  Referring  to 
Pig.  6.1,  the  vector  from  a  general  point  (5;  p;  cp  • )  on  the  trailer 
to  the  point  (x,  r,  0)  can  be  seen  to  be 


D  =  £  (x  -  5'),  r  cos0  -  pcoscp',  r  sin0  -  psin$'^ 


(6.2) 


which  is  the  same  as  (5*2)  except  for  the  primes  on  5 '  and  cp ' . 

A  unit  vector  tangent  to  a  helix  with  an  advance  coefficient  /\(p)  is 


V  p2  +  A2(p)^ 


r*>.  - 


p  simp',  p  costp' 


1 


J 


(6.3) 
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FIG.  6.1  NOTATION  FOR  VELOCITY  AT 

(x,i ;$)  INDUCED  BY  AN  ELEMENT 
OF  A  TRAILER  AT  THE 

TRAILER  ORIGINATES  AT(^^f). 
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She  three  velocity  components  can  be  obtained  by  substituting  (6.2) 
and  (6.3)  in  (6.1) 


00 

u(*,r,e,5,p,f)  -  &  \  (p2  ■  aSBS&tl  -  »)  2L 


2rpoos<p' ”1 


(6.M 


GO 

v(x,r,8,g,p,ep)  .  i-  ]  [(»  -  5')  PPOPT',  -  X(p)>  .lne  tX(p)p  .l^ap 

*+TT  W  |  non  1  <5  /n 


cp '  =cp  [(x  -  5  1  )2  +  r2  +  p2  -  2rp  coscp ' J  ^^2 


00 

w(x,r,  e,5,P,cp)  =  ]— •  J  D^p)r  cosP  - ^lPJP-£0.s2.'-i-(x 

<p'=<p  [(x  -  §')2  +  r2  +  p2  -  2rp  coscp'^J^2 


In  these  equations  the  variable  of  integration  has  been  converted  from  s' 
to  cp'  using  the  relation 


ds'  =  >1  p2  +  X2(p)  dtP'  (6.5) 

The  velocity  components  in  (6.4)  can  be  resolved  into  axial  tangential 
and  radial  components  as  done  before  for  the  sources  and  bound  vortices 


u 


1 


00 

SC  p2  -  pr  cos(<p'  - 

r  n  n 


[(x  -  5 ' )2  +  r2  +  p2  -  2rp  cosfj 


3/2 
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U  *  V  cose  -  V  sine  -  [  U(pLL?.-  fe-l’lp 

cp '=9  [(x  -  5 1  )2+  r2+  p2-  2rpcos  cp'J  3/2 

(6.6) 

m 

^  =  w  Sine  +  V  cose  =  J  L\(p)p  gln(y  -  8>  (x  -  g')p  eosfa'-e)]  &p' 

9 '=9  |jx  -  5’)2+  r2  +  p2  -  2rp  cos9*J  3^2 


The  component  normal  to  the  helicoidal  surface  at  (r,  6)  is  obtained  by 
combining  (6.6)  and  (3.5) 


un(x,r,$,§,p,9)  = 


kn\4Z+  >?(r) 


00 

•  l 


cp'=<p 


rp2-  r2p  cos(9'-e)- X(r)X(p)fr-pcos(9,-e))  -  (x-g*)pX(r)  sin(9'-  e)  dp' 


£(x  -  )2  +  r2  +  p2  -  2rp  cos9'^j 


3/2 


(6.7) 


By  again  introducing  the  approximation 


x  -  5'  =  A(r)e  -  X(p >9'  -  X(p)(9'-  8)  =  -  X(p)u 

which  is  equivalent  to  (3*9)  and  (5*8)  and  replacing  9'-  6  by  v  +  6. 

XL 

(6.7)  becomes 


4tt  /r  +  \2(r) 


i 


£(X(*0X(p)p  -  r2p)  cos(v  +  6k)  +  X(r)X(p)pv  sin(v  +  6fc)  -  r(X(r)X(p)-p2)J  dv 


2,_v.2  .  _2 


2 
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taie  Telocity  Induced  toy  K  unit  trailers  la  obtained  toy 
over  the  blades 


T(r,p,ti) 


a. 

I.S 


t^(x,P,H,k) 


k»l 


, 


(«.9) 


The  total  velocity  Is  obtained  by  Introducing  the  strength  of 
the  trailers  given  in  (4.5),  (4.9)  and  (4.11 )  and  Integrating  over  the 
blade  surface. 


Y  r  _ , 

=  "  J  J  ^  LY^'  **  +  ®)>/p2  +  X2(p)J*(r,P,i»)«i»  dp  + 

rh  0L(p)_8 


(6.10) 


1.0 

|  [v(p,aL(p))  T(r,p,^(p)-e)  “  Y(p,yP))  T(r,p,^j(p)  -  •)  ■  y - j 

rh 

J2  +  >?(P)  dp 


The  first  Integral  In  (6.10)  represents  the  contribution  of  the 
trailers  originating  from  points  in  the  interior  of  the  blades,  while  the 
second  Integral  takes  into  account  the  trailers  starting  at  the  leading 
and  trailing  edges.  The  Cauchy  principal  value  of  the  Integrals  Bust 
toe  taken  whenever  the  range  of  integration  Includes  the  point  (r,6) 
where  the  velocity  is  to  be  determined. 
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While  the  expressions  for  the  relocity  induced  by  the  souroee  end  bound 
vortices  are  not  particularly  sisple,  (6.10)  is  a  inch  acre  difficult 
expression  to  evaluate  due  to  the  fact  that  the  function  T(r,p,u)  involves 
an  integration  over  a  send -infinite  interval.  However,  this  difficulty 
can  be  eliminated  by  considering  T  as  the  sum  of  the  two  functions 


T  =  Tq  +  T ■  (6.11) 

obtained  by  splitting  up  the  range  of  integration  in  (6.3)  in  such  a 
way  that 


TQ(r,p)  =  T(r,p,o) 


w 


>2(r) 


f(v)  dv 


(6.12) 


Tx(r,p,n)  » — \  f(v)  dv  (6.13) 

tor/r2  +  >^2(r)  o 

where  f(v)  is  the  integrand  of  (6.8). 

It  is  evident  that  f(v)  is  an  even  function  if  one  applies  the 
same  argument  used  in  sections  3  wad  5*  Consequently, 

00 

T  (r,») - 1  \  f(»)  4»  (6.1*) 

8n/r2  *  X2(r)  1 

is  one-half  of  the  velocity  Induced  by  a  set  of  X  unit  hellcoldal  vortices 
extending  from  (-*»,•).  Since  the  velocity  field  of  such  a  configuration 
is  two  dimensional  (  it  is  independent  of  s  in  the  (s,n,r)  coordinate 
system)  an  alternate  expression  for  Tq  can  be  obtained  from  a  two 
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dimensional  potential.  A  derivation  of  tfeia  potential  wt  given  bjr 

(a)  \ 

herbs'*'  who  found  the  result  to  be  an  infinite  envies  of  nodi fled 

« 

Bessel  Functions  .  While  this,  in  Itself,  would  attar  no  particular 
computational  advantage  over  the  Integral  r  .■presentation  given  la  (6.13)* 
it  Is  fortunate  that  highly  accurate  asymptotic  approximations  to  the 
sums  of  Bessel  functions  are  available.  The  most  accurate  approxismtlons 
were  developed  by  Wrench^10),  and  it  has  been  found  that  his  formulas  are 
far  more  efficient  than  the  numerical  Integration  required  to  evaluate 
(6.13).  The  exact  and  approximate  expressions  for  Tq,  taken  from  Lerbs ' 
and  Wrench's  work  with  some  modification  in  nomenclature,  appear  In 
Appendix  A. 

We  consider  finally  the  second  part  of  (6. 10),  namely  T^. 

Since  f(v)  is  even,  is  an  odd  function  of  u 

Tx(r,p,-v)  =  -  Tjfopjn)  (6.15) 

This  fact  is  not  only  useful  in  the  numerical  evaluation  of  T^,  but 
also  permits  us  to  draw  some  important  conclusions  regarding  Incidence 
and  camber  which  will  he  considered  in  the  next  section. 


0 

Lerbs’  results  are  expressed  in  terms  of  "Induction  factors” 
which  differ  from  Tq  by  a  constant  factor. 
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7.  Total  Disturbance  Telocity  Bornal  to  the  Blades 

The  total  disturbance  velocity  normal  to  the  first  blade  at 
a  point  (r,9)  Is  obtained  by  sunning  the  effects  of  thickness  (3.17), 
bound  vortices  (5.10)  and  trailers  (6.10) 


un(r,9)  =  u£s)(r,0)  +  u£b\r,6)  +  (7.1) 

We  can  further  subdivide  u^^lnto  parts,  u^°  ^  and  \  where  the 
former  represents  the  contribution  of  Tq  and  the  latter  the  contribution 
of  as  defined  in  (6.12)  and  (6.13).  An  expression  for  u^0^  can  be 
obtained  from  (6.10)  be  replacing  T(r,p,p,)  by  TQ(r,p)  and  bringing  Tq 
outside  the  p,  Integration  (since  it  is  independent  of  n) 


^.(p)-6  _ 

uito)(r)  •  J  T0(^p)j-  J  fy(p,9+e)/p2+  A2(p)]<^  + 

p=rh  L  ‘*=8l(p)_® 

(7.2) 

j^v(aL(p)p)  -  y(^(p)  p)  “—“J  y p2  +  >2(p)|  dp 

This  can  be  sinqplified  using  (4.13)  and  (4.l4)  to 

1.0 

u^o^r)  =  -  J  To(r,p)  dp  (7.3) 

^P 


where  r(p),  as  defined  in  (4.3),  Is  the  total  circulation  around  each 

at  a  radius  p.  This  is  well-known  as  the  "lifting-line"  equation 
which  represents  the  velocity  induced  by  the  trailers  shed  from  K 
concentrated  radial  bound  vortices. 


38 


The  lifting-line  approximation  to  u^(r )  is  obtained  by  deter¬ 
mining  11^°)  from  (7*3)  and-  approximating  u£b)(r,0)  from  two-dimensional 
theory  based  on  the  sections  at  the  particular  radius  in  question.  The 
latter  approximation  is  frequently  referred  to  as  "strip  theory. "  The 
velocity  components  u^S  ^  and  ^  are  assumed  to  be  zero  in  lifting-line 
theory.  However,  these  assumptions  are  valid  only  when  the  aspect  ratio 
of  the  blades  is  large,  which  is  never  true  in  the  case  of  a  marine  pro¬ 
peller.  Consequently,  the  lifting-line  velocity  u^°^(r)  is  of  limited 
usefulness  in  itself  and  should  therefore  be  considered  simply  as  one 
of  the  ingredients  in  the  total  velocity  given  in  (7*1 )• 

We  see  from  (7*3)  that  u^°^(r)  represents  a  disturbance 
velocity  which  is  constant  over  the  chord.  The  remaining  terms  in  (7.1) 
are  generally  functions  of  both  r  and  0  and  depend  on  the  blade  outline 
and  the  load  and  thickness  distribution. 

We  now  consider  the  special  case  where  the  blade  outline, 
chordwise  load  distribution  and  chordwise  thickness  distribution  is 

symmetrical  about  the  lines  0  =  6^  through  the  tips  of  the  blades.  In 

(s) 

this  case,  it  can  be  shown  that  u^  /(r, 9)  is  an  even  function  of  0,  while 

u^(r,0)  and  u^t;L^(r,9)  are  odd  functions  of  0. 

n  11  (s) 

Consider  first  the  velocity  due  to  the  sources,  u^  '(r,e). 

It  was  shown  in  section  3  that  the  function  S(r,p,u)  is  an  odd  function 

of  n*  Recalling  that  u  =  9  -  0,  we  can  write 

S(r,p,cp  -  0)  =  -S(r,p,  -9+8)  (7**0 

Since  the  blade  outline  is  symmetrical,  we  also  know  that 

S^Cp)  -  -  9j(P)  (T.5) 


while  the  fact  that  the  thickness  form  is  synmetrical  about  the  add-chord 
requires  that  the  source  strength  be  an  odd  function  of  cp 


<t(p,<p)  -  -  <y(p-«p)  (7.6) 

The  total  velocity  induced  at  a  point  (r,  6)  by  the  sources 
according  to  (3 .12)  Is 

1.0  _  «t(p> 

P2  +  >2(P)  J*  ®(P><P)  S(r,p,<p-9)  dtp  dp  (7.7) 


«<a>(r,e).  J  / 


while  the  velocity  induced  at  the  corresponding  point  -9  is 


1.0 


Mp) 


u^S\r,-9)  =  J  J 92  +  A2(p)  .J  (r(p,(p)  S(r,p,<p+e)  dtp  dp  (7.8) 

r. 


‘h 

We  next  introduce  the  symmetry  properties  expressed  in  (7*4),  (7-5) 
and  (7.6)  in  (7.8)  to  obtain 


1.0 


■Mp) 


i|ts)(r,-9)  ■  J  Jt?  +  x2(p>  J  »(p,-«p)  S(r,p,-sp-9)  asp  dp  (7.9) 
rh  "^r(p) 

Substituting  -9  for  9  as  the  variable  of  Integration  in  (7 .9)  and 
changing  the  sign  of  the  limits  of  integration  accordingly,  we  obtain 

i-°  _ M>) 


»<S)(r,-8)  =  - 


p2  +  >»2(p)  I  <*(p/p)  S(r,p,cp-e)  dtp  dp  (7.10) 


SfcpCp) 
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Since  this  Is  the  negative  of  (7*7)  except  with  the  upper  and.  lover  Halts 
of  Integration  reversed,  ve  conclude  finally 

-4e)(r,e)  .  «<s)(r,-8)  (T.U) 

33ie  velocities  u^b  ^(r, 6)  and  oan  be  shown  to  be 

odd  functions  of  0  in  exactly  the  same  way.  In  this  case  B(r,p>^,) 
and  T(r,p,p,)  are  odd  functions  of  4  as  was  S(r,p,u)  while  the  strength 
of  the  bound  vortices  and  trailers  is  an  even  function  of  cp.  This 
Introduces  an  additional  minus  sign  in  going  from  the  equivalent  of 
(7*8)  and  (7*9)»  Consequently,  we  conclude  that 


o<b)(r,9)  -  ->4b)(r,-e) 


(T.12) 


u(‘l)(r,8)  -  -  u£bl)(r,-e) 

An  important  conclusion  which  can  be  drawn  from  (7*12)  is  that  a  propeller 
of  zero  thickness  with  symmetrical  blades  and  load  distribution  will  require 
no  additional  incidence  beyond  that  given  by  lifting-line  theory,  i.e., 
the  contribution  of  u^°  \r).  An  additional  incidence  correction  (aside 
from  corrections  due  to  viscosity)  can  appear  only  as  a  result  of  blade 
thickness.  The  lifting  surface  correction  therefore  consists  entirely 
of  camber  in  this  case.  Conversely,  the  thickness  distribution  cannot 
Induce  a  net  camber  since  it  is  an  even  function  of  0. 

These  conclusions  obviously  do  not  hold  if  the  assumed  symmetry 
is  not  present. 
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8.  Hon-Linear  Refinements 


According  to  linear  theory,  the  singularities  representing 
pressure  loading  and  thickness  are  distributed  on  stream  surfaces 
formed  by  the  undisturbed  approach  flow.  Some  improvement  In  accuracy 
would  presumably  be  achieved  if  the  sources  and  vortices  on  the  blades 
were  located  on  the  mean  lines  of  the  sections  at  each  radius  and  if 
the  trailers  were  made  to  follow  the  actual  streamlines  extending  down¬ 
stream  from  the  trailing  edges  of  the  blades.  However,  in  this  case 
the  position  of  the  singularities  would  depend  on  the  disturbance  velocity 
field  which,  in  turn,  would  depend  on  the  position  of  the  singularities. 

While  a  complete  solution  of  this  non-linear  problem  is 
presently  considered  to  be  both  impractical  and  unnecessary,  there  is  one 
refinement  which  may  be  introduced  with  very  little  additional  congplication. 
This  is  customarily  referred  to  as  the  theory  of  "moderately  loaded" 
propellers,  as  contrasted  to  the  strictly  linearized  case  which  is 
termed  "lightly  loaded." 

In  the  moderately  loaded  case,  it  is  still  assumed  that  the 
distortion  of  the  oncoming  flow  due  to  the  radial  component  of  the 
disturbance  velocity  is  negligible.  Consequently,  the  streamlines  remain 
on  cylindrical  surfaces  as  before.  However,  in  this  case,  the  distortion 
of  the  streamlines  due  to  the  axial  and  tangential  disturbance  velocities 
is  taken  into  account  in  an  approximate  way.  It  is  assumed  that  the 
resultant  streamlines  lie  on  helicoidal  surfaces  whose  pitch  includes 
the  lifting  line  disturbance  velocities  obtained  from  (7»3)»  The 
modified  surface  can  be  expressed  in  terms  of  a  hydrodynamic  pitch  angle 
3^,  which  can  be  seen  from  Figure  8.1  to  be: 


(8.1) 


tan  p^r) 


M(*°)(r)  +  yA(r) 
<r  +  t4t°>(r) 


It  should  be  noted  that  the  tangential  disturbance  Telocity 
is  negative  in  Figure  8.1,  which  is  evident  if  one  compares  Figure  8.1 
with  the  velocity  diagram  shown  in  Figure  2.3.  The  plus  sign  in  the 
denominator  of  8.1  is  therefore  consistent  with  this  sign  convention. 

A  hydrodynamic  advance  coefficient  X^  can  be  defined  as 


The  lifting  line  disturbance  velocities  for  a  prescribed  radial 
load  distribution  may  be  obtained  from  (7*3)  as  before.  The  only 
difference  is  that  X  is  replaced  by  in  the  expression  for  Tq.  As  a 
result,  T  becomes  a  function  of  uf^o^  through  A  so  that  (7*3)  represents 
a  non-linear  Integral  equation  rather  than  simply  an  integral.  However, 
a  solution  can  readily  be  obtained  by  an  Iterative  scheme  where  the  J  'th 
approximation  to  u^  is  obtained  by  using  the  value  obtained  in  the 
(j  -  l)bt  iteration  in  determining  Tq.  The  first  approximation  is  simply 
the  linear  one,  i.e.,  -  A. 
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FIG.  8.1  VELOCITY  DIAGRAM  AT  A 
LIFTING  LINE 


Once  \  Is  determined,  the  evaluation  of  the  remaining 
velocity  components  u^s^,  u^5  ^  and  can  he  accomplished  Is  the 

sane  way  as  in  the  completely  linearized  case.  Consequently,  the  only 
additional  complication  Introduced  in  the  theory  of  moderately  loaded 
propellers  is  that  an  iterative  solution  is  required  to  solve  the  lifting 
line  equation. 
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Stannary  of  Expressions  for  Ihdujed  Velocities  at  a  Lifting  Line 

The  axial  and  tangential  velocity  Induced  by  K  helicoidal 
vortex  lines  of  unit  strength  is  given  in  (6.6).  The  lifting-line  case 
is  obtained  by  setting  0  =  cp  =  x  =  0  and  by  replacing  cp'  by  v  + 
and  5  1  by  >(p)v 


aa(r,p)  =  fa 
^(r,p)  =  ^ 


\  {p2  -  pr  cos ( y  +  6Q)  dv _ 

O  |>2(p)v2  +  r2  +  p2  -  2rpcos(v  +  6k)|^2 

(A-l) 

00 

\  (X(p)(r  -  pcos(y  +  6k)  -X(p)pp  sin(v  +  6k)l  dy 

o  [X2(p)v2  +  r2  +  p2  -  2rp  cos(v  +  6fc)J  ^^2 


(9) 

Equivalent  expressions  are  given  by  Lerbsw/  in  terms  of 
modified  Bessel  functions.  With  some  slight  changes  in  nomenclature 
and  sign  conventions,  these  results  are 


for  r  <  p 


un (r, p )  - - - - 

a  4it  X(p) 


l  - 


2KpF^ 

Mp) 


at(r,p) 


■P 


2ttt  A(p) 


where 


nK  )-2& 

A(p) 


(A.2) 
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and  where  JL  and  I  \  are  the  modified  Bessel  functions  of  the  first 


and  second  kind  respectively,  end  the  prime  denotes  differentiation 
with  respect  to  the  arguments. 


For  r  >  p 


ua(r,p) 


5t(r’p) 


-K2p?2 
2tt  X^(p ) 


K 

4ttt 


2KpF, 

"x(7j 


*2  - 


] 

i-K-w1 


nK 


nKp 

X(p) 


n=l 


(A.  3) 


(10) 


The  approximations  to  ?1  and  ?2  developed  by  Wrench'  '  are 


.  .  2  1/4 

i  1  + 

* \  } 

2*^  1  +  q 


2  2 
9<£  +  2  34  -  2 


^ LRj572  oT7)3/2J  w 


^2  ^ 


2Kq  '  1  +  a2' 


2 .1/4 


U-l  24K 


‘  2 

94o  +  2 

(T717)3/2 


3q2  -  2 

^775^ 


log(l 

U-l 


where  -  p/X(p)  4  -  **/XCp) 


(AA) 


(4.(7777-1)  r/ — i  / — gii 

- r —  exp  7i  ♦  <f  -  ✓  1  ♦  ^ 

(  4(7 l  +  ^o2  -  1)  L  J 
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The  quantities  u(r,p)  and  i^(r,p)  can  be  converted  to 
"induction  factors"  as  follows: 


ift(r,p)  =  -  W*-P)  \(r,p)  (A.  5) 

it(r,p)  =  V(r-p)  ut(r,p) 

When  defined  in  this  way,  the  Induction  factors  are  always  positive  for 
a  right-handed  helix.  The  minus  sign  in  the  expression  for  i  (r,p)  is 
not  present  in  Equation  6  of  Reference  (9)  where  induction  factors  are 
first  defined,  but  it  is  incorporated  in  later  expressions  in  the  paper. 

The  reason  for  defining  an  induction  factor  is  to  obtain  a 
quantity  which  approaches  a  finite  limit  as  r-*p.  These  limits,  as 
derived  by  Lerbs  and  others,  are 


ia(r,p)  cosP(p) 
it(r,p)  -  sin0(p) 


(A.6) 


Finally,  the  quantity  TQ(r,p),  which  represents  the  velocity 
normal  to  the  hell  col  dal  surface  at  a  point  r  on  the  lifting  line  can 
be  expressed  in  terms  of  induction  factors  using  (3*5)  and  (A.j) 


TQ(r,p) 


**a(r>P)  -  X(f)  Ut(^P) 

777*5 


r*a(r,p)  ♦  A(r)  it(r,p) 
Mp  -  r)/r2  +  >?(r) 


(A.7) 
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